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Tension dynamics and viscoelasticity of extensible wormlike chains 
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The dynamic response of prestressed semiflexible biopolymers is characterized by the propagation 
and relaxation of tension, which arises due to the near inextensibility of a stiff backbone. It is 
coupled to the dynamics of contour length stored in thermal undulations, but also to the local 
relaxation of elongational strain. We present a systematic theory of tension dynamics for stiff yet 
extensible wormlike chains. Our results show that even moderate prestress gives rise to distinct 
Rouse-like extensibility signatures in the high-frequency viscoelastic response. 

PACS numbers: 61.41. +e, 87.15.La, 87.15. He, 98.75.Da 



Recent experiments have successfully linked the vis- 
coelastic properties of living cells to the rheological be- 
havior of prestressed biopolymer networks p], [1, [E El- 
Single filaments within these networks are well described 
by the wormlike chain (WLC) model, where very stiff 
backbones are idealized as inextensible space curves [f|, 
giving rise to a characteristic / _1 / 2 -divergence of the 
force / required to attain full stretching Q. Here, 
"stretching" needs to be seen as a "straightening" of 
excess length stored in thermal contour undulations Q, 
suggesting the phrase "pulling out stored length" . If sud- 
den forces are applied, stored length can be pulled out 
at first only from growing boundary layers of size £\\(t) 
near the ends due to lon gitu dinal friction with the vis- 
cous solvent @, H, [E [Tol. fll|. The precise time depen- 
dence of £\\ (t) is influenced in a quite subtle way by the 
applied prestress and a non-homogeneous distribu- 
tion of stored length along the contour correspondings to 
a non-uniform tension profile. 

Clearly, as the prestress gets larger, an increasing con- 
tribution to the longitudinal extension stems from the mi- 
croscopic elasticity of the backbone bonds and less from 
thermal undulations. Hence, for an extensible backbone 
the polymer's response is characterized by a local com- 
petition between destroying thermal stored length and 
creating elongational strain. While it has long been rec- 
ognized that stretching modes of long and slender elas- 
tic rods relax extremely fast [l3j |. in a non-equilibrium 
setup they do so only up to the scale where the much 
slower stored length dynamics has equilibrated. Espe- 
cially for bead-spring simulations, where realistically stiff 
backbones often require unfeasibly short time steps, it is 
not immediately clear if and how backbone stretching 
affects the longitudinal relaxation. In this report, we 
present a theory of tension dynamics for stiff yet finitely 
extensible wormlike chains. A brief comparison of exten- 
sible and inextensible polymer models is used to moti- 
vate the ensuing systematic derivation based on the in- 
extensible analog presented in Ref. [8(. We then calcu- 
late viscoelastic response properties and show that even 
moderate prestress can give rise to distinct extensibility 



signatures reminiscent of a Rouse-like dynamics. 

In the WLC model, the polymer backbone is ideal- 
ized as continuous space curve r(s). Contour undulations 
are penalized with a bending energy proportional to the 
squared local curvature. Strict inextensibility would re- 
quire that s be the arclength such that r'(s) is a unit vec- 
tor, and this hard constraint allows exact solutions only 
for special cases In general, Lagrange multipliers of 
varying sophistication are used to enforce miscellaneous 



constraints of different rigidity [li EE EE EE EE El- 
Specifically, for a non-equilibrium scenario with non- 
uniform stored length dynamics, it is inevitable to use 
a local constraint jE IE 11, IE 21 > IE IE > which is in- 
tuitively interpreted as backbone tension. In the case of 
an extensible backbone this tension arises naturally as a 
spring force [24l IE and although it is generally far 



from trivial [27] , our results will permit taking the limit 
from soft to rigid constraints. 

For an extensible but very stiff backbone with only 
small stretching deformations, the Hamiltonian reads 24 1 



n = 



ds[£ p 



k x u 2 



(1) 



where L is the unstretched contour length, £ p is the per- 
sistence length, k x is the stretching elastic constant, and 
u = \r'\ — 1 is the elongational strain. Our theory re- 
lies on the weakly-bending limit r(s) = (s — r\\,r±) T of 
small transverse and longitudinal contour deviations r± 
and ry from a straight line, which gives u ~ —r'^ + \v'J_ to 
leading order. Observing that the polymer's longitudinal 
extension in the limit of large prestress /q 3> k^T j£ v is 
given by 24 1 
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(l-r') = l + (u)-(±r>l), (2) 



we can quantify the simultaneous limits of an only 
slightly extensible backbone and a weakly-bending con- 
tour by requiring that the contributions of longitudinal 
strain (u) = ,fo/ {k-^Tk^) = e x -C 1 and of thermal stored 
length (\r'l) = [k B T£ p /(Af Q )} 1/2 = e th < 1, respec- 
tively, are both small. Although these contributions are 



2 



independent, we will also assume that the force-extension 
Eq. ([2]) is dominated by contour straightening instead of 
backbone stretching, i.e., that (u) -C (^ r '±), which is eas- 
ily fulfilled as long as fo •< / x , where / x = ksTk x ^ 3 /£p^ 3 



is the corresponding crossover force scale [25J, [28[ . This 
assumption is reasonable in most experimental circum- 
stances, considering that f x sw 75 pN [29j for actin, and 
f x « 50 pN for DNA which in fact is close to the 
overstretching transition. For bead-spring simulations, 
however, this condition is much harder to obey [30j |. 
because very small time steps At < A:" 1 are required. 
Also, in special situations like the relaxation from a low- 
temperature initial condition 21| short-time transients 
are quite pronounced, and we will show below that even 
a prestress of only about 0.01/ x gives rise to observable 
effects. 

To quantitatively assess the influence of backbone 
stretching on the dynamics, we proceed with a discussion 
of the equations of motion Cftr = —5H/6r + £, with the 
stochastic noise £ and the free-draining friction matrix 
C = C±[ r ' r ' + C(l — r ' r ')L where £ « 1/2 accounts for 
the anisotropy between transverse and longitudinal fric- 
tion. To leading order in e t h and e x and in the absence 
of external forces, we obtain 

dtr^ = -r"l' + k x (ur' ± )'/£ p + £ ± , (3a) 
(d tr]] + (1 - Qr'xdtrx = -rff" - k x u'/£ p + £y . (3b) 



= I- 1 



Here, we have introduced units such that k^,T 
and C,± — 1, which makes time a length 4 and force a 
length -2 . In the following, we are interested in the pro- 
totypical rheological experiment 0, 0, E| where at time 
t = a small time-dependent force 5f(t) is superim- 
posed on a static prestress fo, which contributes a term 
[fa + Q(t)df(t)][6(L - s) - 6(a)] on the right hand side 
of Eq. (|3bp . In the stationary state at times t < 0, this 
gives k x u/£ p = fo, and the combination k x u/£ p plays the 
role of a tension in Eq. (|3a| also at later times. 

If we assume constant u = £ p fo/k x , we find that the 
transverse part Eq. (|3a[) is correlated on length scales 

*!.(*) With £|~ t 1 /* if t < J-2 and _ (/ Q t)l/2 

if t 3> /q~ 2 @, [HJ. On the other hand, disregarding 
the thermal contribution rii to Eq. (|3b[) for the moment, 
we also find that the diffusive dynamics of the elonga- 
tional strain u is correlated on length scales £ x (t) ~ 
(M/^ P ) 1/2 0. Given now that £ p f Q /k x = C(e x ) < 1, 



it turns out that £ x 3> £j_ except for veryy early times 
t < (£ p /k x ) 2 , where higher order terms become rele- 
vant. Hence, after short initial transients the elonga- 
tional strain u varies slowly with arclength: stretching 
modes relax extremely fast, but only to a local equilib- 
rium value, which is not only nonzero for a prestressed 
filament, but can even show nontrivial large-scale spatial 
variations for the previously mentioned noncquilibrium 
stretching experiments. Thus, elongational strain cannot 
globally equilibrate unless these tension variations have 



propagated through the filament. The latter are linked 
to the dynamics of thermal stored length and therefore 
significantly slowed down by longitudinal friction, and it 
has been shown that the associated characteristic length 
scale ^ii(t) oc e th £x(t) is much larger than the one of 
transverse fluctuations [8( . 

Our goal is now to formulate an equation for the elon- 
gational strain u that integrates over transverse fluctu- 
ations described through Eq. (|3a|) (on the short length 
scale £±), but retains both large-scale spatial variations 
in the tension k x u/£ p (on the scale £») as well as the 
effect of short time transients stemming from the fast 
relaxation of elongational modes. To this end, we em- 
ploy a multiple scale perturbation theory both in space 
and time: small-scale and large-scale spatial coordinates 

l/2_ 

s and s, respectively, account for the different spatial 
correlation lengths of transverse and longitudinal con- 
tour displacements, while slow and fast time variables t 
and t = {eth/s x )t account for long time stored length 
and short time strain dynamics, respectively. Here, the 
condition e x <C £th (resulting from fo <C f x ) is essen- 
tial, and it entails that a r-derivative of u is of order e t h: 
transients from strain relaxation can become comparable 
to thermal stored length. We briefly sketch the analy- 
sis, which in its technical details is quite analogous to 
Ref. [3l|. Taking an s-derivative of Eq. (|3b[) and elim- 
inating rf| = — u + we find to zeroth order in eth- 

—(d t u — d^u—k x d 2 u/£p. Since k x u/£ p = G(fo) while ev- 
erything else is 0(e x ), we find d 2 u — 0, and therefore also 
dgU = dtu = 0. This also implies that Eq. (j3a| becomes 
the linear equation dtr± = —d^r± + k x ud 2 r±/£ p + 
To first order in e t h, Eq. (|3b| now gives (^(d s r±) 2 — 
^e t hd T u/e x = —k x d 2 u/£ p + H(s,s), where H(s,s) sum- 
marizes s-derivatives of terms nonlinear in rj_. These 
vanish upon coarse-graining, i.e., when averaging this 
equation over small-scale variations on the scale £± [3ll |. 
Denoting in such a manner spatially averaged quantities 
with an overbar and replacing r — > t, we obtain the rela- 
tion 



Kd 2 ul 



-(d t {g - u), 



(4) 



where g — ^r' 2 is the local density of contour length 
stored in thermal undulations. This relation, which is 
our main result, formalizes in an intuitive way the op- 
posing effects thermal stored length density g and elon- 
gational strain u have on the backbone tension k x u/£ p : if 
prestress is increased, stored length is destroyed and elon- 
gational strain created, and both lead temporarily to spa- 
tial tension inhomogeneities (curvature), and vice versa 
for decreasing prestress. Further, taking the "inextensi- 
ble" limit e x — » while holding the tension k x u/£n fixed 
leads to the inextensible analog derived in Ref. [8(, and 
the limit Eth — * of a one-dimensional Rouse chain, al- 
though our assumptions cease to hold, nevertheless gives 
a simple diffusion equation for u. 



3 



In order to solve Eq. ((4J, we observe that its bound- 
ary conditions are prescribed through the externally ap- 
plied prestress: u(0,t) = u(L,t) = £ p [f + Q(t)5f(t)]/k K . 
Further, g has to be computed from Eq. (|3ap . which to 
lowest order depends only parametrically on u(s, t). It is 
thus effectively linear in r± and can be solved in Fourier 
space by means of the response function [3l[ x± (Q] t, = 
e - q 2 [ q 2 (t-t')+kAU(t)-u(t'))/i p ]^ where (j^ = Jt d t'u(t') is 

the time-integrated strain. Because spatially averaging g 
over many effectively uncorrelated segments of length £± 
produces an ensemble average, the stored length density 
is given by [^H 



' dq_ 



Solutions to Ec 
Refs. 



xl(g;*,Q) 

q 2 + fo 



+ 2q 2 / dt'xl(q;t,t') 



(5) 



12, 21, 3J. 



^ can now be obtained similarly to 



In the remainder of this paper, we focus on a small 
oscillatory stress Sf(t) = 5 f smut superimposed on a 
largeprestress fo. This situation has been analyzed in 
Ref. [32j for an inextensible filament, and our calculation 
proceeds along these lines. While we cannot actually 
take the linear response limit Sf — > because then the 
assumptions underlyingthe multiple scale perturbation 
theory become invalid [8J, we can still for small enough 
Sf -C fo linearize Eq. (|5|) by writing U(s,t) = £ p fot/kx + 
AU (s, t) and find from Eq. ((U) a simple equation for the 
Laplace transform AU(s, z) [3ll[32T|: 



d§AU(s, z) = M{z)AU(s, z). 



(6) 



The kernel M(z) = C/ 



1/2 



M(zfo 2 )/£ p is defined as 



M(z) = 



'dk 



2k 2 



4fc 4 



fc 2 + l 2fc 2 (fc 2 + l) + z 



A 3 / 2 



(7) 

and (j) x = /o// x *C 1 is the ratio of prestress to the crit- 
ical force / x = k* S /£p 3 . In constrast to the inexten- 
sible case treated in Ref. [j^, where M(z) oc z 1 / 4 for 
all z 3> 1, we obtain now an additional high-frequency 



regime M(z) 



A 3 / 2 



z for z 3> 



Considering that 



M 1 / 2 {z) is the analog in Laplace space to the character- 
istic length scale £\\{t) for the large-scale spatial tension 
variations, we find that £\\(t) ~ {k x t/ '(-p) 1 ^ 2 jlij shows 
a diffusive scaling in the corresponding extensibility- 
dominated short-time regime i < fj, only then crosses 
over to the well-known growth law £\\{t) ~ Ip^t 1 / 8 [Io| . 
before finally arriving at £\\(t) ~ £p 2 f 3 ' 4 t 1 / 2 



11| for 



t^>f z (i.e., I < 1). 



of the crossover time i x = £^J 3 jkx 3 = / x 2 



We emphasize that the scaling 
cannot sim- 
ply be inferred from dimensional analysis, because any 
combination of the lengths £ p and fc^ 1 could be used. 

The observable of main interest is the change SRn (t) in 
projected length Ru , which is through the force-extension 
relation Eq. © related to the integrated change of u — g 



and thus through Eq. ((4]) to the features of the tension 
k x u/ £p\ 



kx 



6R\\(t) = f dt' [ ds(u~g) = -!?-[d s U(L,t)-d s U(0,t)]. 
Jo Jo (£ p 

(8) 

In Laplace space, we thus obtain from a straight- 
forward solution of Eq. © 5R(z) = 3-M 1 / 2 (z) x 

C z 

tanh[|M 1 / 2 (z)]<5/(z), which can be backtransformed in 
the stationary limit t — > oo 32]: 



SRAt) 



J'(w/ - 2 )sin(cjt) - J"(cj/- 2 )cos(^) 



QLf 2 

with the dimcnsionlcss compliances 



(9) 



J'(u>) = -Im 



J"(w) = -Re 



</>c /2 M(iw) tanh J ±</>c /2 M(io>) 



M(\uj) tanh \f \4> c > M{iuj) 



(10) 



Here, <fi c = fo/ fc gives the ratio of fo to the longitudinal 
critical force f c — £ 2 /(( 2 L 4 ) and can be used to distin- 
guish "long" (6 C ^> 1) and "short" (<f) c <C 1) filaments, 



respectively |12|. 

Fig.[T]depicts numerical solutions of Eq. (|10|) for a fixed 
value of <p x = 1CP 2 and different values of </> c , as well as 
results for the corresponding viscoelastic elastic modulus 
(G" + iG") = (J' + U")- 1 . It is straightforward to check 
that Eq. (fTQ|) obeys scaling laws in different intermedi- 
ate asymptotic regimes, which are summarized in Ta- 
bled! Most of these results have been obtained previously 
(cf. Refs. [13, El [H, [H, H HI), and we will therefore 
not comment on them in detail. We merely emphasize 
that for low frequencies to < 1 the presence of a prestress 
leads to the well-known nonlinear response regime with 



its characteristic ^-exponents 



35 



3G] . The intermediate 

regime (?>c'" < ii C 1 with J' ~ J" ~ tD -1 / 2 corre- 
sponds to the previously discussed regime of nonlinear 
tension propagation [8,|Tl| and can only be observed for 
large (j> c 3> 1. Also, for small prestress (<fi c <C 1), there 
is an intermediate regime with J' ~ J" ~ cj" 3 / 4 for 
4>~ 2 » lu » 1 equivalent to the force-free case [22I |33]|. 

The predominant effect of an extensible backbone 
is to produce a Rouse-like scaling in the new high- 
frequency regime 10 f 2 . Here, we find J'(u>) ~ 
J"{uj) - 2 1 / 2 0c /4 0x /4 w~ 1/2 , corresponding to G' ~ 
G" - 2- 3 / 2 ^ 1/4 0x 3/4 w 1/2 . From Fig. □ we conclude 
that a moderate prestress fo of merely 1% of the critical 
force f x suffices to significantly shorten the 7/8- regime 
and to produce a distinct Rouse-like signature, especially 
if cj) c is not too large. 

In summary, we have presented a systematic theory 
of tension dynamics for extensible wormlike chains in- 
cluding the opposing effects of thermal stored length and 



4 



4> c = 100 



- 1 — i — I — i — i — i — I — i — i — i - 




—i i I i i i I i i i_ 



10" 4 10° r , 10 4 

UJ 



L0 



10" 4 10° ,-, 10 4 

UJ 



10° 



FIG. 1: Plot of the moduli J', J" (top) and G',G" (bottom) 
from Eq. JTOj for 4> x = / // x = 10~ 2 and for C = 100 (left) 
and 4> c = 1 (right), where (p c = fo/f c with f c = l 2 /(C, 2 L A ) and 
/x = k 2/3 /£p /3 . Solid lines: J' {&), dashed lines: J" (&'). 
Short lines indicate high-frequency scaling laws of Table U 



(a.) 



J'(u>) 



J"{uj) 



(f>- 2 > uj > 1 
1 > uj > 0~ 1/2 

<?!>c" 1/2 > a 

(b) 



^j 4 uj-y* 
sV'uj- 1 ' 2 



J'(uj) 



UJ > <p~ 

r 2 > uj > <^r 2 
/>c~ 2 > w > 1 

1 > UJ 



,1/4,3/4 „ -1/2 
V X ' W ' 

,1/4 a -7/8 
,1/2 ^ -3/4 



AV2 



TABLE I: Asymptotic scaling results [Trj, [23, [22, [32l, [3J, [23] 
for the compliances J'(a>) and J"(uj) from Eq. (|10[) with t£) = 
tufg 2 for 0c 3> 1 (a) and ^ < 1 (b), respectively, where 
0* = ,fo//x and <j> c = fo/fc with / x = k 2/3 /£p /3 and / c = 
£ 2 /(C 2 L 4 ). 



elongational strain relaxation at short times. These pro- 
duce a Rouse- like scaling in the high-frequency viscoelas- 
tic response, and are expected to be especially relevant 
for the proper design of bead-spring simulations. 
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